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ABSTRACT 


The  development  of  a  speed-of-light  hard-kill  weapon  system  for  military 
applications  represents  a  significant  advancement  in  technology  over  present 
conventional  kinetic  weapon  systems.  Over  the  past  two  decades,  the  US  Navy  has 
successfully  developed  a  megawatt-class  chemical  laser;  however,  under  some  maritime 
environments,  the  high  power  beam  propagation  was  unable  to  delivery  sufficient  energy 
to  kill  a  modem  anti-ship  missile  (ASM)  due  to  significant  atmospheric  absorption  and 
the  resulting  thermal  blooming  process.  A  critical  problem  to  resolve  for  the  shipboard 
high-energy  laser  weapon  systems  is  to  develop  a  shipboard-compatible  megawatt-class 
laser  weapon  at  a  wavelength  where  the  atmospheric  absorption  is  smallest.  The 
megawatt-class  Free  Electron  Laser  (FEL)  has  significant  advantages  over  conventional 
weapon  systems  and  other  chemical  high-energy  laser  systems.  Infinite  magazine,  rapid 
response,  and  wavelength  tunability  make  the  FEL  a  suitable  and  desirable  shipboard 
weapon  system. 

This  thesis  divides  into  four  chapters.  Chapters  I  and  II  introduce  the  FEL  and 
background  theory  of  the  FEL.  Chapter  III  explores  the  analysis  of  the  LANL 
Regenerative  MW  FEL  Amplifier  design  and  optimizes  its  efficiency.  Lastly,  chapter  IV 
summarizes  the  feasibility  of  achieving  the  desired  efficiency. 
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I.  INTRODUCTION 


A.  DEFENSE  AT  THE  SPEED  OF  LIGHT 

In  today's  technology,  weapon  systems  rely  predominantly  on  kinetic  kill,  striking 
a  missile  with  another  missile  or  with  several  bullets.  These  weapons  are  relatively  slow 
and  require  a  longer  detect-to-engage  sequence  when  factoring  in  tracking  time,  launch 
time  and  time  of  flight  to  intercept.  In  searching  for  faster  and  faster  weapons  and  ever- 
decreasing  response  times.  Directed  Energy  Weapons  (DEW)  have  been  on  the  U.S. 
military's  "love  to  have"  list  for  several  decades.  DEWs  are  attractive  for  several  reasons, 
and  the  most  predominant  reason  is  that  they  offer  speed-of-light  defense. 

In  1971,  the  U.S.  Navy  established  the  Navy  High  Energy  Laser  (HEL)  Program 
with  the  intentions  of  developing  speed-of-light  hard-kill  weapon  system.  Initial 
development  focused  on  building  a  CO2  gas  dynamic  laser  technology,  but  a  few  years 
later  its  emphasis  shifted  to  continuous-wave  deuterium  fluoride  (DF)  chemical  lasers 
because  they  propagate  far  better  in  the  maritime  environments.  In  the  early  1980s,  the 
megawatt  class  Mid-Inffared  Advanced  Chemical  Laser  (MIRACL)  was  built  by  TRW 
for  the  Navy  as  the  highest  average  power  continuous  wave  (CW)  laser  in  the  west.  The 
SeaLite  Beam  Director  (SLBD)  was  built  by  Hughes  Aircraft  to  accommodate  and  direct 
the  MIRACL  beam.  The  MIRACL  and  SLBD  have  successfully  tested  at  White  Sands 
Missile  Range,  New  Mexico,  in  engaging  and  destroying  several  types  of  missiles  in 
flight  [1]. 

Although  the  MIRACLE/SLBD  system  was  successfully  tested  against  crossing 
targets  in  a  fleet  defense  scenario,  its  lethality  was  not  sufficient  in  some  ship  self-defense 
engagement  scenarios  due  to  insufficient  crosswind  [1].  Insufficient  crosswind  will  cause 
the  laser  beam  heating  the  air  in  its  own  path  resulting  in  a  defoeusing  optical  lens  known 
as  thermal  blooming  [2].  The  thermal  blooming  effect  can  severely  reduce  the  High 
Energy  Laser  Weapon  Systems'  (HELWS)  lethal  range.  As  a  result,  it's  necessary  to 
replace  the  MIRACL  with  a  megawatt-class  laser  device  at  a  wavelength  where 
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atmospheric  absorption  is  smallest.  An  ideal  solution  is  the  megawatt-class  Free  Electron 
Laser  (FEL)  that  is  tunable  to  the  optimum  propagation  wavelength. 

B.  SHORT  HISTORY  OF  FREE  ELECTRON  LASER 

The  first  generation  of  coherent  radiation  from  free  electron  beams  was 
microwave  tubes  that  relied  on  slow-wave  structures,  and  their  development  received  a 
significant  impact  from  the  radar  development  during  World  War  11.  In  1951,  Motz 
showed  that  an  electron  beam  propagating  through  an  undulator  magnet  can  be  used  to 
amplify  radiation.  After  1960,  research  on  short  wavelength  lasers  has  been  dominant,  but 
they  require  an  atomic  or  molecular  medium  to  operate,  resulting  in  a  limitation  on 
wavelength  tunability.  In  1971,  J.  M.  J.  Madey  of  Stanford  University  brought  attention  to 
what  he  proposed  as  the  'free  electron  laser'.  Relying  on  research  of  synchrotron  radiation 
sources,  Madey  conceived  a  device  that  would  lase  in  the  visible  using  a  beam  of 
relativistic  electrons  as  the  source.  In  1976,  the  first  successful  free  electron  laser  (FEL) 
experiment  was  demonstrated  by  Madey  and  his  co-workers  at  Stanford  University,  in 
that  they  measured  7%  gain  from  an  FEL  configured  as  an  amplifier  at  10pm 
wavelength.  One  year  later,  another  successful  FEL  operation  of  the  same  FEL  configured 
as  an  oscillator.  These  successful  experiments  created  a  large  interest  in  FEL  research. 
Since  then,  scientist  and  experimental  groups  around  the  world  built  FELs  that  operate  at 
frequencies  ranging  fi-om  microwave  to  the  UV  [11].  Three  other  FELs  were  built  soon 
after.  The  first  was  built  at  Laboratoire  pour  I'Utilisation  du  Rayonment 
Electromagnetique  in  Orsay,  France,  which  operated  in  the  visible  range.  The  second 
involved  a  team  from  TRW  at  Stanford  using  the  Superconducting  Accelerator  (SCA), 
lasing  in  the  near  infrared  (IR).  The  third  was  built  at  Los  Alamos  that  lased  in  the  mid- 
infrared  [3].  Now,  there  are  many  FELs  used  for  research  around  the  world. 
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c. 


BASIC  FREE  ELECTRON  LASER  THEORY 


The  PEL  is  a  device  that  transforms  the  kinetic  energy  of  a  relativistic  electron 
beam  into  electromagnetic  (EM)  radiation  [4].  This  is  accomplished  by  sending  a 
relativistic  electron  beam  from  a  linear  accelerator,  microtron  or  a  storage  ring  etc., 
through  an  alternating  magnetic  field  produced  by  a  device  known  as  an  undulator  as 
shown  in  Figure  1.1.  The  interaction  of  the  negatively  charged  electrons  with  the  spatially 
periodic  magnetic  field  induces  a  transverse  acceleration  on  the  electrons  that  causes 
them  to  exchange  energy  with  a  copropagating  radiation  field.  The  kinetic  energy 
extracted  from  the  electron  beam  is  converted  to  electromagnetic  radiation.  The  energy  of 
the  electron  beam  and  the  period  of  the  undulator  magnet  determine  the  wavelength  of 
radiation.  This  system  provides  a  powerful  source  of  tunable,  coherent,  electromagnetic 
radiation  that  can  operate  in  a  wavelength  range  from  about  1  cm  to  the  X-ray  region. 
This  system  can  also  provide  continuous  wave  (CW)  operation,  or  very  short  pulses  in  the 
nanosecond,  picosecond,  or  even  subpicosecond  region. 


Electron  Beam 


Figure  1.1 .  Schematic  diagram  of  a  free  electron  laser. 
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II.  PRINCIPLES  OF  OPERATION  OF  FEL 


A.  ELECTRON  TRAJECTORY  IN  THE  UNDULATOR 


For  better  understanding  of  the  physical  mechanism  underlying  the  operation  of  an 
FEL,  we  first  consider  the  interaction  of  a  free  electron  in  the  undulator  and  an 
electromagnetic  wave.  The  Lorentz  force  equation  and  the  equation  for  the  energy 
exchange  between  the  electron  and  the  electromagnetic  wave  are 


F  =  — =-e| 
dt 


^  -  V 

E+  —  xB 
\  c  J 


d{ymc^) 

dt 


=-ev-E, 


(2.1) 


(2.2) 


where  e,  m  and  v  are  the  charge,  mass  and  velocity  of  the  electron,  respectively;  y  is  the 
relativistic  factor  of  the  electron,  y=l/-^(l-y0^);  fi-v/c;  p  is  the  relativistic 

momentum  of  the  electron,  p='yinv  ;  E  is  the  electric  field  of  the  EM  wave;  B  is  the 
magnetic  field  ;  ymc^  is  the  electron  energy;  and  c  is  the  speed  of  light  in  vacuum  [5]. 

Considering  an  electron  with  velocity  v  =  P  c  and  a  transverse  electromagnetic 
wave  propagating  colinearly,  in  the  z-direction  in  a  vacuum,  then  equations  (2.1)  and 
(2.2)  become 


d(rp) 


dt 


—(E  +  pxB), 
me 


dy  e  -  - 
— P-E. 
dt  me 


(2.3) 


(2.4) 


Assume  the  undulator  has  a  helical  magnetic  field  of  the  form  .B=5[cos(A:„2),sin(A:^z),0], 
where  k^=2n ! is  the  undulator  wavenumber,  and  is  the  undulator  wavelength. 
The  corresponding  circular-polarized,  plane-wave  optical  field  is  then 
E,  =  E  (cos  If/,  -  sin  y/,  0) ,  and  B^ =E(sin  y/,eos  \f/,0 ) , 
where  y/^kz- at-v ^  .  £  is  the  electric  and  magnetic  field  amplitude  in  cgs  units,  t  is 
time,  (j)  is  the  optical  phase,  k=(ol  e  is  the  optical  wavenumber,  and  (O  is  the  optical 
angular  frequency.  Substituting  the  fields  into  Lorentz  force  equations  yields 
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d{YP.)  e 

~dt — 5(-sin(A:^z),cos(A:„z),0)],  (2.5) 
d{yp.)  e 

—  sm^^)  +  B{p^  sm{k„z)- Pycosik^z))],  (2.6) 

and  E(,p^  cosy/- PySini!/),  (2.7) 

where  P  -  Px  +  Pz^  P^  -  Px^  +  Pyj ,  and  p,  =  pS.  For  relativistie  electrons,  p,  «1, 


hence  E{\-  p^)  «  P,B,  so  in  the  right-hand  side  of  (2.5)  the  first  term  is  ignored.  The 
second  term  is 


d(yPJ  eBp^^ 


(-sin(k„z),cos(k^z),0) . 


(2.8) 


dt  me 

Integrating  the  transverse  equation  (2.8)  and  assuming  a  perfect  injection  of  electron 
beam  into  the  undulator  cavity,  which  makes  the  constant  of  integration  zero,  yields 

^  eB^ 


yPi=- 


\mc 


\(cos(k^z),sin(k„z),0) . 


(2.9) 


Divide  both  side  of  (2.9)  by  y  to  get 

Px(^)^  (cos(A:„z),sin(A:„z),0),  (2.10) 

y 

where  the  undulator  parameter  is  K=eB  I  k^mc^  =eB\  1 2mnc^  .  Equation  (2.10) 
describes  the  transverse  motion  of  electrons  in  the  helical  undulator  [6]. 


B.  PEL  RESONANCE  ~  ELECTRON  AND  PHOTON  RACE 

Accelerated  electrons  traveling  at  a  relativistic  speed  along  the  longitudinal  axis  of 
the  undulator  of  the  FEL  radiate  optical  energy  [6].  In  order  for  the  FEL  to  have  gain  (i.e. 
the  fractional  change  in  optical  power  in  a  single  pass  through  the  undulator),  a  net 
transfer  of  energy  from  the  electron  beam  to  the  optical  beam  must  occur.  In  order  for  an 
optimum  exchange  of  energy  from  an  electron  to  the  optical  field  to  occur,  the  electron 
must  oscillate  in  pheise  with  the  optical  field.  This  optimum  energy  exchange  occurs  when 
the  electron  falls  behind  the  optical  fields  by  a  distance  of  one  optical  wavelength  within 
the  span  of  one  undulator  period  [6]. 
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speed  of  electron  ^ 


Speed  of  light  c 


Light  wavelength 


t=0  Undulator  wavelength  X=y  pzC 


Figure  2.1.  Electron-Photon  race. 


Consider  an  optical  wave  that  travels  at  speed  c  and  passes  over  an  electron  traveling  at 
speed  as  shovra  in  Figure  2.1.  At  resonance  condition,  one  wavelength  of  the  optical 
wave  passes  over  an  electron  while  it  travels  through  one  undulator  wavelength  ,  so 

that  light  wavelength  X=  (c  -  p^c)t .  Insert  t=  X^  /  p^c  yielding  X^^—~~Xg .  Substitute 

r  2 

k} 

Pl=  —  into  =  \-  p^  =  \-  P'^  -  Pi  yields 

(\  +  K^)r-^  =\-  Pi .  (2.11) 
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Since  «1,  \-  fi,  » 


\  +  K^ 


1  + 

(2.12) 

2/ 

The  equation  (2.12)  is  the  resonance  condition  that  demonstrates  the  tunability  of  the 
PEL.  By  adjusting  the  undulator  parameter  K,  or  the  energy  of  the  electron  beam 
or  the  undulator  wavelength  the  wavelength  of  the  optical  field  can  be 
tuned  to  a  desired  wavelength  within  a  prescribed  bandwidth  [7].  This  is  a  significant 
advantage  over  conventional  lasers  that  are  typically  constrained  to  a  fixed  wavelength. 


C.  PENDULUM  EQUATION 


This  section  derives  the  equation  of  motion  of  the  electrons  resulting  from  the 
interaction  with  the  optical  field.  Substitute  (2.10)  into  (2.7)  to  get 

eKE 

f  =  —cos(C  +  <P),  (2.13) 

ymc 

where 

^ +  ^  =  y/  +  k^z  =  {k  +  kg)z-(ii^  +  ^.  (2-14) 


For  -7t !  2  <  ^  +  (j>  <  K  !  2  cosine  function  in  (2.13)  is  positive  and  the  electrons  will 
gain  energy  from  the  optical  field,  whereas  for  7t  I2<  +  <!>  <3n  12,  the  cosine  is 

negative  and  the  electrons  will  lose  energy  to  the  optical  field.  Differentiating  (2.1 1)  with 
respect  to  time  t  gives  2y  “^^(l  +  )=  2>9^y0, .  Relating  y  to  p,  gives 


y  _  Y^PzPz 

r 


(2.15) 


Differentiating  (2.14)  once  yields  ^  =^{k  +  kg)z-0)  =  {k-¥k^)cp^-co,  and 
differentiating  twice  yields  ^  =  {k  +  kg)cP , .  Insert  the  resulting  expression  fory?^  into 
(2.15)  to  get 


Y  _ 


(2.16) 
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Since  A  «  then  «  k,  so  (k  +  k^)c»  kc  =  ci).  Using  this  approximation  and  the 
resonance  condition  (2.12),  equation  (2.16)  becomes 


r..  ^ 

r^{l  +  K^)o)  2co/ 


(2.17) 


So, 


f  eKE 

Y  y^mc 


cos(^  +  ^),  then  the  electron  equation  of  motion  is 


f  =— t---‘:os«-  +  «>).  (2.18) 

Y  me 

Equation  (2.18)  describes  the  electron  phase  dynamics  in  the  form  of  the  pendulum 
equation  [6]. 

Consider  time  evolution  through  the  undulator  L/  p^c  for  an  average  electron. 


where  L  is  the  length  of  the  undulator.  Introducing  relevant  dimensionless  time 
T  =  p^ct  /  ct  I L  ,^0500  that  r  will  vary  between  0  and  1  along  the  undulator  length. 


Substituting  dt  =  {L! c')dT  into  the  pendulum  equation  (2.18)  gives 


4'  =  v  =  |a|cos«'  +  (z)),  (2.19) 

o 

where  the  notation  {..)  =  d{..)l  dr ,  \a\  =  AtiNcKLE I y^ is  the  dimensionless  optical 

field  strength,  and  N  =  L/  A^  is  the  total  number  of  undulator  periods.  Equation  (2.19)  is 

the  electron  pendulum  equation  that  governs  the  phase-space  motion  of  the  electrons 
within  the  undulator  and  under  the  influence  of  the  optical  wave  [6]. 
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D.  OPTICAL  WAVE  EQUATION 

An  electromagnetic  wave  will  propagate  in  the  absence  of  a  current  density 
according  to 

1  a'  ^ 


A(r,t)^0, 


(2.20) 


where  the  Laplacian  is  =  a^  +a^  +a^,  r  is  the  vector  position,  t  is  time,  and^  is  the 

optical  vector  potential  [5].  The  optical  electric  and  magnetic  fields  can  be  expressed  in 
terms  of  the  vector  potential  by  the  following  relationships 

£  =  (2.21) 

c  ot 

and  B  =  VxA.  (2.22) 

Assuming  that  the  optical  beam  is  a  circularly  polarized  plane  wave,  then  the  vector 

potential  has  the  form 

-  E(f,t)  ( 1 

"  fsm\i/.cos\u.0i=9le  ^  — 

u 


A  = — ^Y^(sin\|/,cos\[/,0)=9le  \—^(r,t)ee‘^'^ 


(2.23) 


where  ^{r,t)=E(r,t)e'*'^'’'^,  y/ =  kz  -  cot  ■¥  (p  is  the  plane- wave  phase,  e  =  (-/,l,0)  is  the 

polarization  vector,  k  =  2nc  I  co  \s  the  optical  wave  number,  co  is  the  optical  frequency, 
and  <t)(r,t)  is  the  optical  phase  [7].  For  convenience,  we  write 


A  =  —&'(r,t)ee‘°‘ , 
k 

where  a  =  kz-(ot.  Differentiating  (2.24)  with  respect  to  t  gives 

dt  k^  ^ 

where  (..)=d(.)/ dt .  Differentiating  (2.25)  with  respect  to  t  gives 


(2.24) 


(2.25) 


d'^A 


de  k^ 


CO 


(2.26) 


where  (..)=d^(..)/ 6t^ .  Differentiate  (2.24)  with  respect  to  z  gives 
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dz  ^ 


(2.27) 


where  (.)'=d(.)/  dz .  Differentiate  (2.27)  with  respect  to  z  gives 

^  =  +  2ik  &”  -  ,  (2.28) 

where  =d^(..)/  dz^ .  Since  the  laser  beam  consists  of  a  narrow  range  of  frequencies, 
the  variation  in  the  electric  field  amplitude  and  phase  with  respect  to  space  over  an 
optical  wavelength  A  is  small,  so  that  ^'«  and  The  variation  in  the  electric 

field  amplitude  and  phase  with  respect  to  time  over  an  optical  period  is  also  small,  so  that 
^  «6)  S’  and  <fma)<f>  [7].  Therefore,  ^"is  negligible  compare  to  kS' ,  and  S  is 


negligible  compare  to  co  ^ ,  then  (2.26)  can  be  simplified  to 

^  =  ^[-2\oyS-(o^s), 
dt^  k  ^  ’ 

and  (2.28)  can  be  simplified  as 

Substitute  (2.29)  and  (2.30)  into  the  left  hand  side  of  (2.20)  yields 


A(r,t)=  — — 

V\S  +  2ik[s'  +  -s'\ 

dt^  J 

k 

V  c  )_ 

(2.29) 


(2.30) 


(2.31) 


where  Vi  =- 

-L  ^-2  ^,2 


Hence, 


f  ^ 

1 

'^\+2ik 

H - 

v5z 

c  dt ) 

1^=0. 


(2.32) 


The  first  partial  second  derivatives  with  respect  to  x  and  y  describe  how  the  beam  changes 
due  to  diffraction  as  it  propagates.  The  last  partial  derivatives  with  respect  to  z  and  t 
describe  the  propagation  of  a  plane  wave  along  the  z  axis. 

Introduce  dimensionless  time  t  =  ct !  L  so  that  r  varies  from  0  to  1  along  the 
undulator  length.  It's  convenient  to  follow  the  wavefront  with  a  new  coordinate 
z  =z-ct .  With  this  coordinate  change. 
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A-  _ 

dz  \dz  j 


dz  \  d  ( ^  ^ 

dz  \dz j  dr  8z' 


(2.33) 


since  dz  /dz  and  dr/ dz  =  0 .  Also, 


ii-  =  i 

c  dt  c 


dz^ 

\dt  j 


d  ( dx^  d 


dz 


\dt  J 


dr 


A  lA 

dz  L  dx’ 


(2.34) 


since  dz  /dt  =  -c  and  dx / dt  =  c / L .  Therefore,  the  propagation  operator  in  (2.32) 
becomes 


dz  c  dt  dz 

Hence,  the  full  wave  operator  is 


did  d  i 

— +  — +  1 - + 

dz  L  dx  J 


d  1  5^ 

L  dx 


(2.35) 


V{+2ik\ 


ill- 


d  1  d' 
_+ - 

c  dt 


+2ik 


lA 

Ldx 


(2.36) 


Introducing  the  dimensionless  transverse  variables  3c  =  xylk/2L  and  y  =  y^kl2L  ,  so 


that  -^  =  {kl2L)^  and  ^  =  {k/2L)^.  Hence  V 
OX  ox  oy  ^ 


=  {k/2L{ 


d 


2  \ 


\dx^ 


_  ^  q2  q2 

{k / 2Z)V^ ,  where  = — y  + ' — t  ■  Substituting  into  (2.36)  and  (2.32)  yields 


Sc" 


=  0. 


(2.37) 


4  "  dx_ 

Equation  (2.36)  is  known  as  the  parabolic  wave  equation  [7].  In  dimensionless 


parameters,  the  optical  field  is  represented  by  a  =  a  e'^ 


where 


|a|  =  Att  NeKLE  /(y  )  qc  £ ,  so  ^  and  a  differ  by  a  multiplicative  constant.  Therefore, 

exchange  for  a,  the  dimensionless  parabolic  wave  equation  (231)  becomes 


a  =  0. 


(2.38) 


4  dx_ 

In  the  presence  of  the  eleetron  beam  current  density,  the  optical  wave  equation 
(2.20)  is 
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v72  1  ^  4;r- 

V  — Ar,0  = - 


(2.39) 


where  is  the  transverse  electron  beam  current  density  [8].  Substituting  (2.30)  into 


(2.39)  yields 


...f  ^  1 

—  V,  +2ik  —  + - 

k  [dzcdt 


- Ji- 


(2.40) 


Performing  a  dot  product  with  e*  =  (/,1,0)  on  both  sides  of  (2.40)  and  using 
e*  -  6  =  ■  (-i,l,0)  =  2 ,  yields 


2  \dzcdtj 


(2.41) 


The  total  transverse  electron  beam  current  is  the  sum  of  all  single  particle  currents.  The 
single  particle  current  for  one  electron  in  the  beam  is  =  -ecP^^5^^\x-r-)  [8]  where 

is  the  position  of  the  electron  and  is  the  three-dimensional  Dirac  delta- 


function.  The  total  beam  current  is 


(2.42) 


From  (2.10),  the  transverse  motion  is  K/ y\cosk^z,sink^z,(!))- 

^He  {(-  K  / y)iee~'^°‘ }.  For  convenience,  we  write  ,  so  that 

Y 

-i(k„z+a) 

=iecK(e'  ■e)y- - S^^\x-r,) .  (2.43) 

/  r 


We  can  simplify  the  sum  ^  by  assuming  a  local  electron  density  p,  so  that  (2.41) 


becomes 


dz  c  dt 


1  2  -ll  ^  ^  5P  -)  •  ^ 

-V,+ik—\ — —  &  =  -iTTieKpi  — 


(2.44) 
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where  ^  =  k^z  + a  is  the  electron  phase  in  the  combined  undulator  and  optical  fields, 
and  (...)  is  a  weighted  average  of  sampled  electrons  at  a  specific  site  which  replaces  the 

summation  over  all  electrons  within  one  optical  wavelength.  Define  the  dimensionless 
beam  current  j  =  SN(e?rKLy  pl(y^mc^) ,  the  wave  equation  becomes 

^^1  +  ;—  a(x,y,z,T)=-j(e-‘^).  (2.45) 

In  some  cases  where  the  optical  wave  diffraction  is  not  significant,  we  can  ignore  in 

(2.45),  then 

a  =  -j{e-‘^).  (2.46) 

The  optical  wave  equation  (2.46)  couples  the  electrons  to  the  copropagating  optical  wave 
by  relating  the  electron  phase  to  the  change  in  the  optical  field.  The  equation  (2.46)  can 
be  expressed  in  terms  of  amplitude  and  phase, 

=  -y(cos(4'  +  ^)} ,  (2.47) 

and  i>  =  /.  (sin(4'  +  (^)) .  (2.48) 

rl 

Equation  (2.47)  reveals  that  when  the  electrons  become  bunched  near  (4'  +  ^)  « ;r,  the 
optical  field  amplitude  will  grow,  and  when  bunched  near  (4'  +  ^)  »  0 ,  the  optical  field 
amplitude  will  decrease.  Equation  (2.48)  describes  the  optical  phase  evolution. 
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E.  GAIN  DEGRADATION  DUE  TO  ELECTRON  BEAM  QUALITY 


In  the  design  of  an  PEL,  there  is  often  a  trade-off  between  high  beam  current  and 
good  beam  quality  as  determined  by  the  accelerator.  The  dimensionless  current 
j  X  IN'^2}''^  where  I  is  the  actual  beam  current  ,  N  is  the  number  of  undulator  periods, 
and  X  is  the  optical  wavelength.  Usually  the  optical  wavelength  A  is  fixed  for  a  given 
application  so  that  j  is  maximized  by  adjusting  current  /  and  the  undulator  periods  N. 
Increasing  I  tends  to  degrade  the  beam  quality  from  the  accelerator,  whereas  increasing  N 
tends  to  increase  the  PEL  sensitivity  to  the  beam  quality  because  of  narrowing  the  gain 
spectrum  bandwidth  [9]. 

Emittance  s  describes  one  aspect  of  electron  beam  quality.  A  real  electron  beam 
has  an  emittance  s  describing  how  electrons  enter  the  undulator  with  a  range  of  angles 
6^  with  respect  to  the  undulator  axis,  and  at  a  range  of  transverse  positions  slightly 
off  axis.  The  electron  beam  may  also  have  a  small  energy  spread  txyly .  The  energy 
spread  and  emittance  of  the  electron  beam  translate  into  a  spread  in  phase  velocities  in 
phase-space.  The  phase  \  clocity  spread  due  to  energy  spread  is  given  by  [10] 

/Xv  =AkN  is.y  !  y ,  (2.49) 

and  the  phase  velocity  spread  of  an  electron  beam  due  to  emittance  is  given  by 

(2.50) 

1  +  A 

Electrons  enter  the  undulator  with  a  spread  in  phase  velocities  due  to  a  spread  in 
energies,  and  are  randomly  distributed  in  both  the  longitudinal  and  transverse  directions. 
The  typical  distribution  of  electron  phase  velocities  is  assumed  to  have  a  Gaussian  form 
about  initial  phase  velocity  , 

/<,(9)  =-/=—.  (2.51) 

where  -AnN  b^y  ! y  is  the  standard  deviation  and  q  is  the  phase  velocity  away  from 
Vg  distributed  [9].  The  characteristic  function  F^jfr)  of  the  distribution  /c(q)is  defined 
as 
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W  =  ffo(i)e-‘”dq. 

Evaluating  the  above  integral  equation  for  a  Gaussian  yields 

If  the  beam  quality  is  perfect  so  that  cr^  =  0 ,  then  (r)  =1 .  However,  when  the  beam  is 
not  perfect,  the  characteristic  function  |F^(r')[  decays  over  the  length  of  the  undulator  in 
a  characteristic  time  on  the  order  of  l/cr^  [9]. 

Poor  beam  quality  reduces  the  single-pass  PEL  gain  by  degrading  electron 
bunching.  As  an  example,  consider  an  PEL  design  with  N  =  25  periods  having  a 
dimensionless  current  density  of  j  =10  and  initial  phase  velocity  of  =3  operating  with 
an  initial  dimensionless  optical  field  Og  =  1.  In  Pigure  2.2,  electron  beam  quality  is 
assumed  to  be  nearly  perfect,  i.e.  cr^j  =1 .  Some  electrons  gain  energy  and  move  ahead  of 

the  average  flow,  while  other  electrons  lose  energy  to  the  radiation  field  and  move  back 
behind  the  average  flow  causing  visible  spatial  bunching. 


Pigure  2.2.  Phase  space  evolution  for  a  nearly  perfect  electron  beam  quality. 
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In  Figure  2.3,  the  longitudinal  energy  spread  increases  to  cTq  =6,  and  this 
degradation  in  beam  quality  causes  much  less  bunching  when  comparing  to  Figure  2.2. 
The  poor  beam  quality  is  responsible  for  a  decrease  in  final  gain  to  G/=  0.10534  and  a 
smaller  optical  phase  shift  ^ . 


Figure  2.3.  Phase  space  evolution  for  low  gain  with  electron  beam  degradation. 

F.  TAPERED  UNDULATOR 

The  characteristics  of  the  FEL  interaction  can  be  altered  by  changing  the 
undulator  properties,  i.e.  undulator  polarization,  wavelength,  or  field  strength,  along  its 
length  [5].  As  an  FEL  reaches  saturation,  the  electrons  have  given  up  enough  energy  so 
that  they  are  no  longer  in  resonance  with  the  optical  wave,  which  decreases  the  gain. 
Tapering  changes  the  undulator's  resonance  condition  along  its  length  so  that  the 
electrons  remain  in  resonance  longer,  hence  increasing  the  overall  gain  in  strong  fields 
[7].  From  the  resonance  condition  (2.12),  in  order  to  extract  energy  out  of  the  electron 
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beam,  i.e.  reducing  y,  and  convert  into  radiation  field  and  still  maintain  resonance 
condition,  we  can  either  taper  the  undulator  wavelength  1^,  i.e.  reducing  X^,  or  taper  the 

undulator  parameter  K,  i.e.  reducing  the  undulator  field  strength  B.  However,  tapering  K 
is  easier  since  the  xmdulator  gap  needs  to  be  increased  to  decrease  B.  Tapering  essentially 
has  the  same  mathematical  result  as  accelerating  the  electrons.  This  artificial  acceleration 
is  represented  by  a  dimensionless  parameter  6.  A  tapered  undulator  represents  an 
advantage  over  the  untapered  undulator  design  when  the  artificial  acceleration  6  exceeds 
the  deceleration  that  can  be  obtained  without  taper  in  strong  optical  fields.  When  the 
undulator  wavelength  is  decreased,  the  artificial  acceleration  d  is  given  by 
5 »  -IttNAX^  /  X^ .  When  the  undulator  field  strength  is  decreased, 

S  -4;iNK^AB/ B{\  + K^) .  With  tapering  is  included  in  the  undulator,  the  pendulum 
and  wave  equations  become 

oo  o 

^  =  V  =  5 +  \a\  cos(^  +  <!>) ,  (2.52) 

O 

and  a  =  -j  <  e~‘^  >, 

where  ^k^(T')dT' +  lcz [10].  The  wave  equation  remains  the  same,  but  the 

simple  pendulum  equation  includes  an  additional  constant  torque  due  to  the  artificial 
acceleration  S . 

I  I  1  |l/2 

The  necessary  and  desirable  criteria  for  tapering  is  |«|  >  ^  -^rl  - 

condition  on  the  left  is  for  trapping  electrons  in  the  tapered  phase  space;  the  middle 
condition  is  that  the  tapered  acceleration  exceeds  the  natural  untapered  deceleration;  and 
the  condition  on  the  right  is  that  the  tapered  PEL  works  in  the  strong  field  regime  [10]. 

Electrons  initially  near  the  phase  4^^  «0  are  accelerated  by  both  the  strong  optical 

field  |a|  and  the  artificial  "torque".  The  electrons  initially  near  undergo  a 

cancellation  between  the  optical  field  and  the  torque,  leaving  them  trapped  in  closed 
orbits.  The  taper  is  effective  because  electrons  near  the  phase  for  gain  Q  ’^n  trapped, 
and  continue  to  contribute  to  the  interaction.  Electrons  near  the  phase  for  absorption 
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^  »0  are  taken  away  from  resonance  and  eventually  stop  interacting.  The  imbalance 
leads  to  a  net  gain  [10]. 

Compeired  to  the  untapered  undulator,  the  tapered  undulator  is  more  efficient  in 
strong  optical  fields,  but  has  smaller  gain  in  weak  optical  fields.  The  efficiency  of  a 
tapered  undulator  is  estimated  at  rjg  =5l(%nN).  For  example,  if  ^«280;r  for  an 

undulator  with  23%  taper  in  wavelength  over  N  =  300  periods,  then  the  efficiency  is 
estimated  at  rjg  k1Q/(2N)  which  increases  by  a  factor  of  70  over  the  natural  efficiency 

t]  «l/(2iV)  [10].  Tapering  is  a  design  to  obtain  higher  overall  gain  and  efficiency  in 
strong  optical  fields. 

An  example  of  a  tapered  undulator  in  strong  optical  fields  is  shown  in  Figure  2.4 
with  an  untapered  undulator  shown  in  Figure  2.5.  Both  designs  are  based  on  an  undulator 
length  of  Z  =  6  m  corresponding  to  A^=  300  periods.  The  dimensionless  current  density  is 

j  =  12000  »  K  in  the  high  gain  regime,  and  initial  field  =150  » ;r  in  the  strong  optical 

field  regime.  In  Figure  2.4,  the  tapering  rate  S  =  260;r  turns  on  at  =0.19  along  the 
undulator  resulting  in  a  final  gain  of  Gf  =  556  and  FEL  extraction  efficiency  of 
r|«12.5%.  Initially  the  gain  increases  exponentially  until  =0.19  where  the  taper 

turns  on  and  saturation  is  reached.  With  tapering,  the  gain  continues  to  increase  slowly 
instead  of  decreasing  after  saturation,  as  shown  in  the  upper  right-hand  portion  of  the 
Figure  2.5.  The  phase  velocity  of  the  untrapped  electrons  near  v  *740  is  shifted  by  the 
tapering  acceleration  S  whereas  the  trapped  electrons  remain  near  resonance  due  to 
strong  optical  fields.  The  lower  right-hand  portion  of  the  Figure  2.4  is  the  optical  phase 
shift,  and  this  phase  shift  w  1 .4;r  over  the  length  of  the  undulator  is  a  characteristic  of 
high  gain. 


19 


Figure  2.4.  Phase  space  evolution  for  high  gain  in  strong  fields  with  a  tapered  undulator. 


Figure  2.5.  Phase  space  evolution  for  an  untapered  undulator. 
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III.  SIMULATIONS  OF  THE  REGENERATIVE 
MW  FEL  AMPLIFIER 

A.  MW-CLASS  FEL  CALCULATIONS 

In  order  to  destroy  an  incoming  missile,  the  laser  power  required  on  the  missile 
surface  is  about  10  kW/cm  over  a  5  cm  radius  spot  for  a  two  second  duration.  The 

extinction  coefficient  due  to  aerosols  at  sea  level  is  «  «  0.03  km‘\  and  through  e~‘^ 
describes  the  absorption  loss  at  X.»lpm  wavelength  over  a  distance  z.  Let  P,  be  the 
Transmitted  Power  that  must  leave  the  ship  to  destroy  a  missile  at  a  range  It,  then  the 
power  at  the  target  is  Pj=P,e~'^.  If  the  range  of  the  missile  is  5  km,  the  laser  power 

which  must  leave  the  ship  is  then  P,  =  [10  kW/cm^]  [7r(5  cm)^]  »  910  kW. 

As  an  example,  consider  an  FEL  with  N  =  300  periods  using  a  beam  of  100  MeV 
electrons  with  a  peak  current  of  400  A,  then  the  peak  electron  beam  power  P  =  400 
Ax  100  MV  «  40  GW.  The  micropulses  generated  by  the  FEL  have  a  length  of  0.001  m 
and  are  separated  by  6  m,  then  the  fraction  D  =  0.001/6  denotes  the  duty  cycle  of  the  FEL. 

The  average  power  generated  within  the  FEL  is  given  by  P=PxD«6.7  MW.  The 
efficiency  necessary  to  destroy  an  ASM  is  ti=P,  /  P«13%.  Therefore,  the  single-pass 
FEL  efficiency  r)  needed  to  supply  the  Laser  Power  to  destroy  the  incoming  missile  is 
13%.  An  approximation  of  the  natural  single-pass  efficiency  of  an  FEL  is  only 
t]  =  1/2N  =  2%  so  that  a  tapered  undulator  with  high  efficiency  is  needed.  The  tapered 
undulator  described  in  the  previous  section  would  be  more  than  adequate.  This  type  of 
illustration  is  the  essential  driving  factor  in  determining  the  FEL  parameters  for  possible 
MW-class  weapon  systems. 
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B.  LANL  REGENERATIVE  MW  EEL  AMPLIFIER 


The  US  Navy  is  investigating  the  possibility  of  using  a  MW-class  Free  Electron 
Laser  for  ship  self-defense  against  anti-ship  missiles  (ASM).  A  design  workshop  resulted 
in  two  possible  MW  FELs,  the  oscillator  and  the  regenerative  amplifier.  The  MW 
oscillator  design,  although  more  compact  overall,  presents  a  challenge  in  recirculating  an 
intense  electron  beam  because  of  a  phenomena  known  as  coherent  synchrotron  radiation 
(CSR)  feedback  [11].  The  regenerative  amplifier  design  uses  a  longer  undulator  and 
relies  on  higher  extraction  efficiency  to  achieve  high  average  power  so  that  less  current  is 
needed  from  the  accelerator  and  no  recirculation  is  required.  The  regenerative  amplifier 
design  has  been  proposed  by  LANL  and  Boeing  with  the  analysis  being  performed  at 
NFS.  This  section  explores  the  analysis  of  regenerative  amplifier  design  and  optimizes  its 
efficiency. 

A  simple  schematic  of  the  regenerative  amplifier  design  is  shown  in  Figure  3.1. 

Fee€l)ackLoq) 


Electron  Dunq) 

Figure  3.1.  Conceptual  MW  class  Amplifier  FEL. 


The  amplifier  design  does  not  use  energy  recovery  so  there  is  no  requirement  to 
bend  the  high  energy  beam  which  can  cause  deterioration  of  the  electron  beam  quality. 
However,  without  energy  recovery,  higher  extraction  efficiency  is  required.  One¬ 
dimensional  computer  simulations  are  used  to  describe  a  single  pass  of  the  optical  wave 
and  to  optimize  the  efficiency  using  a  tapered  undulator.  The  tapered  xmdulator  represents 


22 


an  advantage  when  the  electron  phase  acceleration  S  exceeds  the  deceleration  that  can  be 
obtained  without  taper  in  strong  optical  fields.  The  single-pass  PEL  efficiency  7  of  a 
tapered  undulator  is  the  fraction  of  the  electron  beam  energy  converted  to  laser  light.  The 
limit  on  rj  is  determined  by  the  maximum  taper  rate  6  that  can  maintain  trapped 
electrons  over  the  number  of  periods  N  in  the  undulator  [6].  Gain  G  is  the  fractional 
change  in  power  of  the  optical  field  in  a  single  pass  through  the  undulator.  In  weak  fields 

without  tapering,  the  PEL  amplifier  has  high  gain  described  by  G{t)  « 

[12].  The  gain  is  exponential  in  t  along  the  undulator  with  growth  rate  proportional  to 
where  j  =  %Nie7t KLf pF ! y^mc^  is  the  dimensionless  current  density, 
/>=3  X 1 0’  I  {A)  /  ecjtr^  is  the  electron  beam  density,  and  F  =  "area  of  electron  beam'V'area 
of  light  beam"  is  the  filling  factor. 

Table  3.1  shows  the  proposed  parameters  of  the  MW  regenerative  amplifier 
developed  at  LANL. 


Parameters 

Proposed  Value 

Beam  energy  E 

100  MeV 

Beam  radius 

0.17  mm 

Energy  spread  Ay/y 

0.02% 

Pulse  duration 

20  ps 

Pulse  frequency 

500  MHz 

A 

Peak  current  I 

400  A 

Average  current  1 

0.2  A 

Undulator  parameter  K 

1.71 

Undulator  length  L 

6  m 

Number  of  imdulator  periods  N 

300 

Undulator  period 

2.0  cm 

Optical  wavelength  X 

1.0  pm 

Table  3.1.  Parameters  for  the  MW  regenerative  amplifier  design  developed  at  LANL 
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1.  ID  Efficiency  Optimization 

The  proposed  MW  RAFEL  uses  a  beam  of  100  MeV  electrons  with  a  peak  current 
of  400  A,  yielding  a  peak  electron  beam  power  of  40  GW.  The  average  current  of  the 
RAFEL  is  0.2  A,  so  that  the  average  electron  beam  power  is  20  MW.  The  single-pass 
FEL  efficiency  77  is  the  fraction  of  power  extracted  from  electron  beam  in  one  pass 
through  the  undulator.  The  RAFEL  will  feedback  on  the  order  of  0.01%  to  1%  of  the 
optical  power  and  require  an  extraction  efficiency  of  approximately  10%  to  15%  to 
provide  2  MW  to  3  MW  of  optical  power  in  the  infrared  (IR). 

In  an  attempt  to  achieve  the  desired  efficiency,  the  initial  undulator  length  is  I  = 
4m  corresponding  to  N  =  200  periods.  The  dimensionless  initial  optical  field  and 

tapered  undulator  of  strength  S  are  varied  to  find  the  optimum  efficiency.  The  maximum 
extraction  efficiency  found  was  8%  which  is  less  than  the  desired  efficiency. 

In  an  effort  to  improve  the  design  efficiency,  the  length  of  the  undulator  is 
increased  from  4  m  to  6  m.  The  change  in  undulator  length  corresponds  to  dimensionless 
current  density  of  j  « 12,000  with  N  =  300  undulator  periods.  The  values  of  and  S 
are  varied  to  find  the  optimum  efficiency  using  numerical  simulations.  The  electron 
beam  has  a  Gaussian  spread  =1.0,  and  the  field  strength  is  =150  with  taper  rate 

S  =260^  starting  at  time  =0.19  along  the  undulator.  Figure  3.2  is  the  3D  plot  of 
efficiency  versus  and  S  .  The  maximum  efficiency  77^3,^  is  13.3%  with  gain  G  =  556. 
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The  effect  of  changing  away  from  its  optimum  value  is  smaller  than  changing  5  away 
from  its  optimum  value  as  shown  in  Figure  3.3. 


Efficiency  Optimization 


For  a  better  understanding  of  the  physical  mechanism  imderlying  the  operation  of 
this  MW  RAFEL,  examine  Figure  3.3  closely  in  five  regions.  Region  I  is  the  "optimum 
region".  As  shown  in  Figure  3.4,  the  taper  turns  on  at  -  0.19,  and  approximately  half 

of  electrons  were  trapped  near  resonance  =  0 ,  resulting  in  an  acceptable  efficiency  of 

7=13.3%.  The  final  gain  Gy- =556  corresponds  to  feedback  of 

/=l/(Gy.+l)=1.8xl0“\  The  large  optical  phase  shift  (j)  in  the  lower  right-hand  portion 

of  Figure  3.4  represents  good  optical  guiding  which  is  an  important  part  of  a  high  gain, 
high  current  FEE. 


Figure  3.4.  Phase  space  evolution  for  the  MW  RAFEL  at  optimum  efficiency. 
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Region  II  is  the  "high  tapering  rate  region".  Increasing  the  tapering  rate  S  beyond 
the  optimum  value,  the  efficiency  drops  to  8.3%  due  to  fewer  electrons  being  trapped  near 
resonance  =  0  as  shown  in  Figure  3.5.  The  taper  turns  on  at  =  0.19  as  before.  The 
final  gain  Gy  =341  corresponds  to  feedback  of  /=l/(Gy^ +l)=2.9xl0■^  ^  is 
relatively  small  which  would  cause  less  optical  guiding. 


-n/2  t  3n/2  Of  1 


Figure  3.5.  Phase  space  evolution  for  the  MW  RAFEL  with  high  tapering  rate. 
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Region  III  is  the  "low  tapering  rate  region".  For  S  below  the  optimum  value,  the 
efficiency  drops  to  6.5%  due  to  insufficient  taper  rate  as  shown  in  Figure  3.6.  While  many 
electrons  are  trapped  near  resonance  =  0 ,  the  low  taper  rate  does  not  extract 
significant  energy  from  them.  The  taper  again  turns  on  at  =0.19.  The  final  gain 
Gf  =  294  corresponding  to  feedback  of  / = 1  /(Gj-  + 1) = 3 .4  x  1 0“^ 


Figure  3.6.  Phase  space  evolution  for  the  MW  RAFEL  with  low  tapering  rate. 


Region  IV  is  the  "high  optical  field  region".  Starting  with  high  optical  field 
strength  ,  the  separatrix  does  not  have  a  chance  to  grow  and  trap  electrons,  resulting  in 
a  lower  efficiency  of  9.5%  as  shown  in  Figure  3.7.  The  taper  turns  on  at  =  0,16 .  The 
final  gain  =144  corresponds  to  feedback  of  / = 1  /(Gj^  +1) =6. 9x10"^. 


Figure  3.7.  Phase  space  evolution  for  the  MW  RAFEL  with  high  optical  field  strength. 


30 


Region  V  is  the  "low  optical  field  region".  Starting  with  low  optical  field  strength 
Qg ,  the  initial  separatrix  is  too  small  so  that  more  electrons  leak  out  from  the  separatrix, 
resulting  in  a  lower  efficiency  of  10.9%  as  shown  in  Figure  3.8.  The  taper  turns  on  at 
=  0.32 .  The  final  gain  Gf  =  25,700  corresponds  to  feedback  of 
/=1/(G^+I)=3.9xl0-'. 


Figure  3.8.  Phase  space  evolution  for  the  MW  RAFEL  with  low  optical  field  strength. 


The  above  results  for  efficiency  should  be  considered  an  over-estimate  since  the  effect  of 
diffraction  of  the  optical  beam  or  longitudinal  pulse  effects  are  not  included. 
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2.  3D  Diffraction  Optimization 

From  Chapter  II,  the  parabolic  wave  equation  can  be  written  as 

a(x,y,z,T)=-j{e-'^).  (3.1) 

Equation  (3.1)  governs  the  dynamics  of  the  optical  wave  over  many  optical  wavelengths  in 
the  longitudinal  z-dimension,  and  the  operator  properly  describes  the  diffraction  of  the 
optical  wave  in  the  transverse  (x,>')  directions.  A  coherent,  freely-propagating  optical 
wave  of  radius  Wo  will  spread  due  to  natural  diffraction  as  it  traverses  the  undulator  length 
L.  The  Rayleigh  length,  nwlIX,  is  the  characteristic  distance  over  which  the  initial 

optical  mode  area  doubles  in  size.  If  the  optical  mode  area  becomes  much  larger  than  the 
electron  beam,  the  optical  field  will  decouple  from  the  electron  beam.  Therefore,  it  is 
important  that  the  Rayleigh  length  is  comparable  to  the  undulator  length  L,  and  the  optical 
mode  waist  Wo  is  not  too  much  larger  than  the  electron  beam  radius  rj,  [13]. 

Three-dimensional  simulations  are  used  to  describe  a  single  pass  of  the  optical 
wave  through  the  undulator.  Initially,  the  optimal  values  from  one-dimensional  simulation 
were  used  as  the  inputs  for  a  3D  simulation.  The  resulting  3D  simulation  is  shown  in 
Figure  3.9  with  =  1 50 ,  S  =  2607r ,  y  =  2.6  x  1  O'* ,  and  the  taper  turning  on  at  =  0.2 . 

The  evolution  of  the  optical  mode,  ja(x:,r)|  shown  in  the  upper-left  window  of  Figure  3.9, 
is  dominated  by  diffraction  after  the  taper  is  turned  on  at  =0.2.  The  intense  electron 

beam  provides  guiding  over  the  first  one  fifth  of  the  undulator,  but  after  tapering  begins 
diffraction  overcomes  this  focusing  over  the  rest  of  the  undulator.  The  top-center  window 
shows  a  cross-sectional  view  of  the  optical  mode,  |a(x,3^)|,  at  the  end  of  the  undulator, 

and  this  optical  mode  has  a  drop  in  optical  power  in  the  center.  This  drop  is  due  to  a  loss 
of  guiding  of  the  optical  mode  [14].  The  window  in  the  top-right  lists  the  dimensionless 
parameters  for  this  simulation;  cr^  and  correspond  to  the  radial  size  of  the  electron 

beam,  while  cr^  and  cr^  account  for  beam  quality. 
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FEL  WAVEFROKTS 


Figure  3.9.  3D  diSraction  simulation. 

The  graph  at  the  right-middle  shows  the  development  of  the  optical  phase  along  the 
undulator  length.  At  =  0.2 ,  the  optical  phase  change  along  the  undulator  length 

becomes  negative  and  drops  rapidly  which  causes  the  optical  wavefront  to  be  excluded 
from  the  electron  beam  decreasing  the  interaction  strength.  The  lower-right  window  is  the 
plot  of  the  dimensionless  optical  power  growth  along  the  undulator,  and  the  plot  of  the 
natural  logarithm  of  the  single-pass  gain  along  the  undulator.  At  the  middle-left  shows  the 
bunching  current  in  the  electron  beam  development  along  the  undulator,  with  the  end  view 
in  the  center.  The  bottom-left  plot  is  the  electron  beam  phase  velocity  evolution  along  the 
undulator,  and  at  the  bottom  center  shows  the  electron  phase  space  plot  at  the  end  of  the 
undulator.  Initially  the  optical  guiding  confines  the  optical  mode  until  the  taper  turns  on  at 
=  0.2 ,  then  diffraction  effects  begin  to  spread  the  beam.  The  optical  phase  begins  to 

grow  until  the  taper  turns  on,  then  the  change  in  optical  phase  is  negative  as  the  guiding 
effect  is  lost.  The  bottom-center  window  shows  that  at  the  end  of  undulator,  no  electrons 
are  trapped  at  resonance,  and  the  tapering  effect  has  failed. 
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Figures  3.10  through  3.14  show  the  dependence  of  efficiency 77  as  and 

5  are  varied.  These  two  parameters  can  have  a  significant  impact  on  efficiency 
when  moved  away  fi'om  their  optimum  values.  The  two  graphs  of  each  figure  are 
plotted  with  the  same  data  with  the  top  one  is  in  3D  and  the  bottom  one  is  in  2D. 


3D  Efficiency  Optimization  (for  a  ==100) 


Figure  3. 1 1.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  tum-on  time  at  ao  =  100. 


Increasing  optical  field  strength  Oo  beyond  its  optimum  value,  the  separatrix  does 
not  have  a  chance  to  grow  and  trap  electrons,  resulting  in  a  slightly  lower  efficiency  of 
7. 1%  as  shown  in  Figure  3.14. 


3D  Efficiency  Opfinnization  (for  a„=250) 


Contcni  r  Plot 


Figure  3.14.  3D  and  Contour  plots  of  Efficiency  vs.  Tapered  Field  Strength 
and  Taper  tum-on  time  at  ao  =  250. 
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Next,  explore  the  effect  on  efficiency  by  varying  ao  and  at  a  fixed  value  of  ^ . 
As  shown  in  Figure  3.15,  starting  at  a  low  value  of  the  taper  field  rate  5  =  S57c ,  the 
maximum  efficiency  obtained  was  3.2%  with  «0.2.  This  low  efficiency  is  due  to 
insufficient  taper  rate. 


3D  Efficiency  Optimixation  (for  5=55tc) 


Figure  3.15.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  ^  =  55;^ . 
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As  d  increases,  the  efficiency  increases  to  at  a  slightly  lower  t,  as 

shown  m  Figure  3,16.  The  peak  efficiency  increases  to  7^^  «7.2%  corresponding  to 
optimum  values  of  ^  =  180;r,  ao  =225,  and  an  earlier  taper  tum-on  time  of  =  0. 18 . 


3D  Efficiency  Optimization  (for  S^lOSit) 


Figure  3. 16.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  turn-on  time  at  <5  =  105;t 


As  shown  in  Figure  3.17,  when  the  optimum  value  of  taper  rate  is  used,  8  -  180;r , 
varying  the  field  Uo  has  little  effect  on  efficiency,  whereas  changing  the  taper  turn-on  time 
from  its  optimum  value,  » 18 ,  results  in  a  significant  efficiency  drop. 


3D  Efficiency  Optimiziition  (for 


Contour  Plot 


Figure  3.17.  3D  and  Contour  plots  of  Efficiency  vs.  Optical  Field  Amplitude 
and  Taper  tum-on  time  at  ^  =  180;r . 
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Now  consider  the  three-dimensional  simulation  of  a  single  pass  of  the  optical  wave 
through  the  undulator  at  the  optimum  values  of  ^  =  180;r,  =  0,18,  and  ao  =  225.  As 

shown  in  Figure  3.19,  optical  guiding  occurs  and  stays  with  the  intense  electron  beam 
through  the  undulator  as  seen  in  upper-left  and  middle-left  windows.  The  top-center 
window  shows  a  much  tighter  cross-sectional  end  view  of  the  focused  optical  power  than 
seen  in  Figure  3.9.  The  bottom-center  window  shows  approximately  50%  of  the  electrons 
are  trapped,  resulting  in  optimized  gain  and  efficiency.  The  middle-right  window  shows 
the  optical  phase  initially  increasing  until  the  taper  turns  on  at  «  0.18 ,  then  decreasing 
slowly.  The  FEL  performance  in  Figure  3.19  is  much  improved  over  that  of  Figure  3.9.  In 
steady-state  operation,  the  final  gain  0^=616  corresponds  to  feedback  of 

/  =  l/(G^+l)  =  1.6xlO“^ 


rsL  WAVsntoHis 


=  0. 18  Ti  =  0.072  8  =  1807C  Gf=616 

Figure  3.19.  3D  diffraction  simulation  at  optimum  values 
S  -  ISOtt ,  =  0. 18 ,  and  Oo  =  225. 
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IV.  CONCLUSIONS 


The  development  of  a  speed-of-light  hard-kill  weapon  system  for  military 
applications  may  represent  a  significant  advancement  in  technology  over  present 
conventional  kinetic  weapon  systems.  The  PEL  has  the  potential  for  use  as  a  weapon 
system  against  anti-ship  and  theatre  ballistic  missiles.  The  infinite  magazine,  rapid 
response,  and  wavelength  tunability  at  high  powers  make  it  suitable  and  desirable  for 
shipboard  self-defense  and  for  land-based  theater  ballistic  missile  defense. 

Three  dimensional  computer  simulations  indicate  that  the  proposed  LANL 
Regenerative  MW  PEL  design  does  not  yield  the  desired  extracted  efficiency  Q^  t]k  15% . 

The  optical  field  amplitude  ao,  taper  rate  5 ,  and  taper  tum-on  time  were  varied  to 
search  for  the  optimum  efficiency,  resulting  in  maximum  efficiency  of  only  rj^  «  7.2% . 
Changing  5  or  away  from  their  optimum  values  has  a  significant  effects  on  the  MW 

RAPEL  design  efficiency,  but  changes  in  Oo  do  not  have  much  of  an  effect.  The  efficiencies 
found  here  are  still  an  over-estimate  since  they  do  not  include  longitudinal  pulse  effects.  It 
remains  to  be  determined  whether  a  smaller  efficiency  is  acceptable  or  if  15%  efficiency  is 
actually  necessary  in  the  PEL  weapon  design.  If  15%  efficiency  is  necessary,  a  longer 
undulator,  perhaps  8  or  10  meters  may  be  sufficient  to  reach  that  goal. 
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